Tuning the Josephson current in carbon nanotubes with the Kondo effect 
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We investigate the Josephson current in a single wall carbon nanotube connected to supercon- 
ducting electrodes. We focus on the parameter regime in which transport is dominated by Kondo 
physics. A sizeable supercurrent is observed for odd number of electrons on the nanotube when the 
Kondo temperature Tk is sufficiently large compared to the superconducting gap. On the other 
hand when, in the center of the Kondo ridge, Tk is slightly smaller than the superconducting gap, 
the supercurrent is found to be extremely sensitive to the gate voltage Vbg- Whereas it is largely 
suppressed at the center of the ridge, it shows a sharp increase at a finite value of Vbg- This increase 
can be attributed to a doublet-singlet transition of the spin state of the nanotube island leading to 
a TT shift in the current phase relation. This transition is very sensitive to the asymmetry of the 
contacts and is in good agreement with theoretical predictions. 
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Metallic single wall carbon nanotubes have attracted a 
lot of interest as ID quantum wires combining a low car- 
rier density and a high mobility. Depending on the trans- 
parency of the interface between the nanotube and the 
electrode, the conduction ranges from insulating behav- 
ior and strong Coulomb blockade at low transparency to 
nearly ballistic transport with conductance close to Ae^ /h 
when the transparency is high^ ^. The intermediate con- 
duction regime is particularly interesting, because in the 
case of an odd number of electrons on the nanotube a 
strongly correlated Kondo resonant state can form, where 
the magnetic moment of the unpaired spin is screened by 
the spins of the electrons in the leads^. Moreover, when 
the carbon nanotubes are in good contact with supercon- 
ducting electrodes, it is possible to induce superconduc- 
tivity and observe supercurrents, as was first investigated 
in ungated suspended devices^ ^ . Proximity induced su- 
perconductivity was then explored in gated devices with 
evidence of a strong modulation of subgap conductance, 
but in most cases no supercurrent was observed^ ^ ^ ^ . 
More recently, tunable supercurrents could be detected 
in the resonant tunneling conduction regime, where the 
transmission of the contacts approaches unit}'^^^^. In 
this regime, the discrete spectrum of the nanotube is 
still preserved and the maximum value of supercurrent 
is observed when the Fermi energy of the electrodes is at 
resonance with "the electron in a box states" of the nan- 
otube. A superconducting interference device was fabri- 
cated with carbon nanotubes as weak links: supercurrent 
TT phase shifts occurred when the number of electrons in 
the nanotube dot was changed from odd to even^^ corre- 
sponding to the transition from a magnetic to a nonmag- 
netic state. Sharp discontinuities in the critical current at 
this — TT transition in relation with the even-odd occupa- 
tion number of the nanotube quantum dot were also ob- 
served in single nanotube junction device^. As pointed 
out in the SQUID experiments^ a Josephson current can 




FIG. 1: a) Scanning electron micrograph of the electrical on- 
chip environment of the nanotube. b) Cross section of the 
nanotube with electrodes and backgate c) Schematic of the 
carbon nanotube Josephson junction with its environment ac- 
cording to the extended RCSJ model. 



be observed in the Kondo regime, when the Kondo tem- 
perature Tk is large compared to the su perconducting 
gap A, confirming theoretical prediction j-'^^ l -'^^ l -'^^ l -'^^ l -'^^ l -'^^ l 
and previous experiment j^ ^ ^^ ^ (with no determination of 
supercurrents though). 

In the present work we explore in detail this competi- 
tion between Josephson and Kondo physics by monitor- 
ing on the same device as a function of the gate voltage 
the bias dependence of the differential conductance in 
the normal state and the Josephson current in the super- 
conducting state. The value of this current is precisely 
determined by fitting the data with a theoretical model 
explicitly including the effect of the electromagnetic envi- 
ronment onto the junction. From the normal state data, 
we extract the charging energy U and the Kondo temper- 
ature Tk as a function of the gate voltage, and calculate 
the sum of the couplings to the electrodes F = F^ + F^^ 
from the known form of T^^^. The asymmetry Vr/Vl 
is determined from the value of the conductance on the 
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Kondo resonance for T <C Tk- The value of this cur- 
rent is precisely determined by fitting the data with a 
theoretical model explicitly including the effect of the 
dissipative electromagnetic environment onto the junc- 
tion. As expected a supercurrent is observed (for an odd 
number of electrons on the nanotube) when Tk > A 
and when the asymm etry of the tran smission of the elec- 
trodes is not too largd^Smmillllll^ On the other hand, 
when Tk < A, the magnetic spin remains unscreened at 
all temperatures leading to a 7r-j unction with a very low 
transmission of Cooper pair^^. We have particularly ex- 
plored in this paper the intermediate regime where the 
Kondo temperature can be tuned with the gate voltage, 
on a single Kondo ridge, from a value slightly below the 
superconducting gap at half filling to a value Tk > A. 
A sharp increase of the supercurrent is then observed, 
which is related to the transition from a magnetic dou- 
blet state to a non-magnetic singlet state of the nanotube 
island. We compare these results with functional renor- 
malization group (FRG) calculations for the single impu- 
rity Anderson model^^ in a wide region of gate voltage. 

We grow single wall carbon nanotubes (SWCNT) by 
chemical vapor deposition on thermally oxidized, highly 
doped silicon wafers. Individual SWCNTs are located 
relative to predefined markers and contacted to Ti/Al 
leads using electron-beam lithography (Fig. [l^,b). The 
Ti/Al leads are superconducting below IK. Measure- 
ments were done in a dilution refrigerator with a base 
temperature of T = 40mK. The cryostat was equipped 
with a three stage filtering system consisting of LC filters 
at room temperature, resistive micro-coax cables and fi- 
nally micro-ceramic capacitors in a shielded metallic box, 
which also contained the samples and was tightly screwed 
onto the cold finger. As additional filters against voltage 
fluctuations, resistors of the order of 2kQ were imple- 
mented on-chip as meanders in the Au lines connect- 
ing the superconducting electrodes to the contact pads 
(Fig. [l^). Measurements were done with a lock-in tech- 
nique with either a AC- volt age of 5 to 10 /iV, measur- 
ing differential conductance dl/dV^ or a AC-current of 
lOpA, measuring differential resistance dV/dl^ both as 
a function of an additional DC-bias voltage or current, 
respectively. Fig. [T]^ shows a schematic view of our de- 
vice. It consists of a Josephson junction with current 
/j(0) parallel to a capacitor Cj and a shunt resistance 
Rj. We obtain a rough estimate of Cj ~ lOOaF from 
the charging energy U = /2Cj. The junction resis- 
tor Rj represents the contribution of quasi-particles to 
dissipation at frequencies of the order of the plasma fre- 
quency of the junction uop. The outer capacitance C 
represents the capacitances of the leads to ground within 
the electromagnetic horizon of the junction, which we 
estimate as ^-Kcjuj^ to some centimeters, depending on 
Rj. It is therefore mainly determined by the capaci- 
tance of the metallic leads to the highly doped backgate, 
and can be estimated to C 8pP^. The series resis- 
tances R represent the lithographically defined on-chip 
resistors. This situation corresponds to the "extended 
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FIG. 2: Color scale plot showing gate and voltage bias de- 
pendence of the differential conductance, a) normal state, b) 
superconducting state. White lines correspond to vertical sec- 
tions of G versus V^-ias 

in the middle of the four Kondo ridges 
labeled A to D. c) gate dependence of the Kondo tempera- 
ture extracted on the Kondo ridges far from the degeneracy 
points. Red lines are fits to eq. [l] the black line corresponds 
to Tc — l.OiT. IJ and F are given in units of A. 



resist ively and capacitively shunted junction" (extended 
RCSJ) model^o 2^ which we wih use later. 

We first characterize the CNT-quantum dot with 
the electrodes driven normal by a small magnetic field 
(lOOmT). The color scale plot of dl /dV as a function of 
Vuas and Vbg (the "charge stability diagram"), displays 
a regular sequence of "Coulomb blockade diamonds" over 
a wide range of Vbg- In other regions this pattern is re- 
placed by a smoother gate dependence characteristic of 
Fabry-Perot oscillations^. For further analysis, we con- 
centrate on the gate voltage region between 3 and 4 Volts 
(see Fig. |2|i). It shows a fourfold periodicity in the size 
of the Coulomb blockade diamonds^ which indicates a 
clean nanotube with the twofold orbital degeneracy of the 
electronic states preserved. We extract U = 2.5±0.3meV 
as estimated from the size of Coulomb blockade diamonds 
with an odd number of electrons. In all states with odd 
occupation, the Kondo effect manifests itself through a 
high conductance region around zero bias, the so called 
Kondo ridge. The Kondo temperature Tk can be esti- 
mated from the half width of the peaks of these lines 
which can be fitted by Lorentzian curves^^. The Kondo 
temperature goes through a minimum of the order of IK 
in the center of the ridges and increases on the edges. 
It is possible to follow this gate dependence along the 
Kondo ridges A to C (Fig. The intensity on ridge 

D is too weak for such an analysis. Tk can be well fitted 
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FIG. 3: a) Color scale plot of the differential resistance as 
a function of DC bias current and backgate voltage, b) I-V 
characteristics for different backgate voltages. The red lines 
are fits using equation^ c) Conductances Gj — 1/Rj (black 
circles), Gs (blue line) and Gn (red line), d) Gate voltage 
dependence of the critical current extracted from fits to eq. 

m 



by the expression predicted by the Bethe Ansatz^B 



Tk = ^/UT/2exp 



8UT 



|4e^ - 



(1) 



where e is the energy shift measured from the center 
of the Kondo ridge. Taking /7=2.5meV as determined 
above, the value of Tk at e = leads to the characteris- 
tic coupling energies F = Fi^ + r^ between the electrodes 
for each Kondo ridge (Fig. The gate voltage Vbg 

dependence of Tk yields the ratio a between the electro- 
static energy cVbg and the Fermi energy of the nanotube 
e equal to 20 ± 1. This value agrees to within 20% with 
the value deduced from the normal state conductance 
data. This fit based on the single level Anderson impu- 
rity model is only valid when U/T is sufficiently large 
and |e| <C U. Since for all peaks T <C T^, the maximum 
conductance of the ridges yields the asymmetry of the 
coupling Tr/Tl, 6.8, 6.2, 2.5 and 70 for A, B, C and D, 
respectively. 

By switching off the magnetic field, we allow the leads 
to become superconducting. Fig. shows dl /V for the 
same gate voltage range as in Fig. |2^ but for a smaller 
bias voltage range. The BCS-like density of states in the 
electrodes leads to new features in the stability diagram, 
which are horizontal lines at Vbias = i2Ao and ±Ao 
due to the onset of quasiparticle tunneling and Andreev 
reflection, respectively. We can derive the value of the 
superconducting gap Aq = 0.15 ± 0.02meV, which cor- 
responds well to the expected Tc — IK for the bilayer 
Ti/Al. Although they have similar Kondo-temperatures 
in the normal state, the Kondo-ridges A, B, C and D show 
a very different behavior when the electrodes are super- 
conducting. The Kondo ridges A, B and D are suppressed 



or reduced in amplitude by superconductivity, reflecting 
that Tk < Tc in the center of the Kondo ridge. Contrary 
to ridge A,B and D, that show a minimum at zero bias, 
there is a strong enhancement of conductance in ridge C, 
with G reaching a value roughly four times larger than in 
the normal state. Note that the states A, C and D show 
an enhancement of conductance at Vbias = A/e, to val- 
ues larger than the conductance at Vbias = 2A/e. This 
enhancement of the first Andreev process is due to the 
"even-odd" effect in A ndreev transport, that has recently 
been described irP^. 

For a measurement of the supercurrent the device has 
to be current biased. We simultaneously use AC and DC 
bias while measuring the resulting voltage drop across 
it. From the AC part, we obtain data on the differen- 
tial resistance (Fig. |3^). By numerical integration we 
get I-V curves that show a supercurrent branch and a 
smooth transition to a resistive branch with higher resis- 
tance (Fig. jsj)). The transition between the two regimes 
is not hysteretic, and the supercurrent part exhibits a 
non zero resistance Rs at low bias even if we subtract 
the contribution of the on-chip resistances R. This be- 
havior is common in mesoscopic Josephson junctions that 
have a high normal state resistance of the order of the re- 
sistance quantum h/e'^. To extract the supercurrent, we 
use a theory that explicitly includes the effect of the dis- 
sipative electromagnetic environment onto the junction 
in th e frame of the already mentioned extended RCSJ- 
modefl^Il^^. Using the external resistor R'^^ (Fig. ^) 
and temperature T as input parameters, we can thus ex- 
tract the critical current Ic and the junction resistance 
for every measured backgate voltage, from a fit to : 
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where r] = hVbias/'^eRkBT and laix) is the modified 
Bessel function of complex order cP^. The resulting val- 
ues are plotted as a function of Vbg in Figs. (Sj) and d. 
The junction conductance Gj = 1/Rj relates well to the 
differential conductance Gs = l/Rs extracted from the 
AC part of the current biased data, especially around 
the resonance degeneracy points where the conductance 
is high. The normal state conductance Gn deviates from 
Gs most notably in the states B and C, where Kondo 
physics plays a key role. In state B, Gn > Gs whereas 
the ridge in state C persists in the superconducting state, 
resulting in a further enhancement of the conductance. 

The supercurrent exhibits peaks at the maximum val- 
ues of the conductance (Fig. [sji). Its behavior between 
peaks varies strongly in states with even and odd occu- 
pation. In the following we focus on the Kondo ridges 
B and C. Whereas the supercurrent on ridge C varies 
nearly proportionally to the conductance, we observe a 
sharp drop of supercurrent on ridge B as illustrated on 
Fig. [sji (arrows). This can be understood consider- 
ing that Tk is smaller on ridge B than on ridge C (see 
Fig. [2]^). As a result the magnetic moment of the ex- 
cess electron on the nanotube remains unscreened in the 
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FIG. 4: Comparison with FRG calculations: Calculated /c(e) 
in units of Iq = eA/h for F = 2A, F/U = 0.11, Fr/Fl = 6 in 
a) and F = 2A,r/[/ = 0.2, Fr/Fl = 3 in b). Experimental 
data for /c(e) are shown for Kondo ridge B in a) and Kondo 
ridge C in b). 



superconducting state near e = which is not the case 
for ridge C. At this stage we compare our experimental 
findings to approximate zero-temperature FRG calcula- 
tions which allow for extracting both the stability regions 
of the screened singlet and magnetic unscreened dou- 
blet phases and the complete supercurrent-phase relation 
/j(^) within a model of an Anderson impurity coupled 
to two superconducting electrodes^^. It was previously 
observed that the Josephson current cannot be simply 
described as a single function of the ratio T^/A^^. The 
relevant parameters are the on site Coulomb repulsion 
energy t/, the level position e related to Vbg^ and the 
transmission of the electrodes F^^ and F^, compared to 
the superconducting gap. In Fig. [4] we present a com- 
parison between our experimental data to the theoreti- 
cal predictions for Ic (defined as the maximum value of 
over as a function of the position in energy of 
the Anderson impurity level, for similar values of theoret- 
ical and experimental parameters. It is remarkable that 
the gate dependence of the critical current on the Kondo 



ridges can be qualitatively described. That is in partic- 
ular, the existence of a singlet, doublet (O/tt) transition 
on ridge B at nearly the right value of gate voltage when 
renormalized with the charging energy U (arrows), and 
the absence of a (O/tt) transition for ridge C correspond- 
ing to higher values of Tk- One point of disagreement 
between theory and experiment concerns the amplitude 
of the critical current in the doublet state (tt junction) 
region (around Vbg = 3.65V). Ic is theoretically found 
to be reduced by only a factor of two compared to its 
value in the singlet region whereas basically no trace of 
superconductivity could be detected experimentally. It is 
known however, that the 7r-phase current computed from 
the approximate FRG is too large compared to numerical 
renormalization group (NRG) data, which are known to 
be more accurate but only available at the center of the 
Kondo ridg^I^. 

Finally, let us emphasize the importance of the asym- 
metry of the transmission of the electrodes which tends 
to reduce considerably the supercurrent. This is par- 
ticularly striking for the data on ridge D, for which 
^r/^l — 70, where no supercurrent could be measured 
in spite a value of Tk of the order of IK. This reduction 
of Ic by the asymmetry of contacts is also found in FRG 
calculations^^. It is moreover accompanied by a modifi- 
cation of the stability regions for the singlet (screened) 
and doublet (magnetic) states which strongly depend on 
the phase difference between the superconducting elec- 
trodes. The non-magnetic singlet state is stabilized with 
respect to the magnetic doublet state in the vicinity of 
= TT, which results in a strong modification of the 
current phase relation which unfortunately cannot be 
checked in a single Ic measurement^. 

In conclusion, we have shown that it is possible to tune 
the amplitude of the supercurrent in a carbon nanotube 
Josephson junction within a Kondo ridge in a very narrow 
range of gate voltage. Our data are in good qualitative 
agreement with theoretical findings and should stimulate 
new experiments where the whole current phase relation 
is measured^^. We acknowledge fruitful discussions with 
Y. Avishai, S. Gueron, W. Belzig, A. Levy-Yeyati, T. 
Novotny, J. Paaske, and B. Rothlisberger. This work 
was supported by the EU-STREP program HYSWITCH 
and by the Deutsche Forschungsgemeinschaft via FOR 
723 (CK, VM). 
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